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ABSTRACT 

Context. M 82 is one of the nearest and brightest starhurst galaxies. It has been extensively studied in the last decade and hy now it is 
considered as the prototypical extragalactic PDR and a reference for the study of the star formation feedback. 

Aims. Our aim is to characterize the molecular chemistry in M 82 at spatial scales of giant molecular clouds (GMCs), ~100 pc, to 
investigate the feedback effects of the star formation activity. 

Methods. We present interferometric observations of the CN 1—»0 (113.491 GHz), N 2 H+ 1—>0 (93.173 GHz), H(41)a (92.034 GHz), 
CH 3 CN (91.987 GHz), CS 3^2 (146.969 GHz), C-C 3 H 2 3 i, 2 ^ 22 ,i (145.089 GHz), H 2 CO 2 o, 2 ^ 1 o,i (145.603 GHz) and HC 3 N 16-»15 
(145.601 GHz) lines carried out with the IRAM Plateau de Bure Interferometer (PdBI). PDR chemical modelling is used to interpret 
these observations. 

Results. Our results show that the abundances of N 2 H^, CS and H'^CO^ remain quite constant across the galaxy confirming that 
these species are excellent tracers of the dense molecular gas. On the contrary, the abundance of CN increases by a factor of ~3 in 
the inner x2 bar orbits. The [CN]/[N 2 H^] ratio is well correlated with the H(41)a emission at all spatial scales down to ~100 pc. 
Chemical modelling shows that the variations in the [CN]/[N 2 H+] ratio can be explained as the consequence of differences in the local 
intestellar UV field and in the average cloud sizes within the nucleus of the galaxy. 

Conclusions. Our high-spatial resolution imaging of the starburst galaxy M 82 shows that the star formation activity has a strong 
impact on the chemistry of the molecular gas. In particular, the entire nucleus behaves as a giant photon-dominated region (PDR) 
whose chemistry is determined by the local UV flux. The detection of N 2 H+ shows the existence of a population of clouds with 
A,,>20 mag all across the galaxy plane. These clouds constitute the molecular gas reservoir for the formation of new stars and, 
although distributed all along the nucleus, the highest concentration occurs in the outer xl bar orbits (R~280 pc). 

Key words. Galaxies: Individual: Messier Number: M 82, Galaxies: Nuclei, Galaxies: Starburst, ISM: Molecules, Molecular 
Processes, Radio Lines: Galaxies 


1. Introduction 

M 82 is one of the nearest and brightest starburst galaxies. 
Located at a distance of 3.9 Mpc, and with a luminosity of 
3.7xl0'° Lq, it has been the subject of continuum and line ob¬ 
servations at all wavelengths from X-rays to the radio domains. 
Several molecular line studies indicate that the strong UV field 
has heavily influenced the physical conditions, kinematics and 
chemistry of the interstellar gas (Mao et al. 2000; Weiss et al. 
2001a,b; Garcia-Burillo et al. 2001, 2002; Fuente et al. 2005, 
2006, 2008; Aladro et al. 2011). In fact M 82 has become the 
prototypical extragalactic PDR and a reference for the study and 
interpretation of star formation feedback in extreme starbursts 
near and far. 

The presence of a stellar bar in M 82 has been invoked as a 
mechanism to fuel the star formation activity in the inner r ~ 30" 


Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany), and IGN (Spain). 


(500 pc) disk of the galaxy. The stellar bar, studied in detail by 
Greve et al. (2002), could have formed during the encounter of 
M 81 and M 82. The orbits sustaining the bar potential, denoted 
as Xl orbits, are oriented parallel to the bar major axis, and ex¬ 
tend out to r ~ 30" (500 pc) in the disk. The M 82’s stellar 
bar seems to have room for X 2 orbits, which are oriented paral¬ 
lel to the bar minor axis, in the inner r ~ 5"(~90 pc). Evidence 
supporting the existence of X 2 orbits in M 82 is derived from 
the position-velocity diagrams obtained in several tracers of the 
ISM, including CO, HI and [Nell] emission line data (Shen & 
Lo 1995; Neininger et al. 1998; Wills et al 2000; Achtermann & 
Lacy 1995). Coupled stellar population synthesis and photoion¬ 
ization models reveal two major episodes of star formation over 
the past 10 Myr (Forster-Schreiber et al. 2003). The first episode 
(~ 10 Myr) took place throughout the central regions of M 82 and 
was particularly intense at the nucleus while the second episode 
(~5 Myr) occurred predominantly in a circumnuclear ring at a 
radius of ~90 pc and along the stellar bar. This inside-outside 
scenario can be understood as resulting from the gravitational 


1 




D. Ginard et al.: The chemical footprint of the star formation feedback in M 82 on scales of ~100 pc** 


Table 1. Observations 


Line 

Freq. 

(GHz) 


beam 

(") 

Fields 


Tft/Flux 

K/(Jy/beam) 

Date or ref ^ 

C‘»0 

1^0 

109.782 

3.8' 

'x3.5" 

PA-4° 

( 0 ", 0 ") 


7.6 

WOl (Zero-spacing) 

CO 

2^1 

230.538 

1.5' 

'xl.4" 

PA0° 

(-30", - 8 "), (-20", 

- 6 ") 

11 

WOl (Zero-spacing) 







(-10", -4"), (10", 4") 









(20", 8 "), (30", 10") 



CN 

1^0 

113.491 

2.5' 

'x2.3" 

PA 170° 

( 0 ", 0 ") 


17 

May-Dec 2005 

H'^CO* 

l -»0 

86.754 

5.9' 

'x5.6" 

PA 105° 

( 0 ", 0 ") 


5 

GOl 

HCO 

1^0 

86.670 

5.9' 

'x5.6" 

PA 105° 

( 0 ", 0 ") 


5 

GOl 

HOC* 

1^0 

89.487 

4.4' 

'x3.6" 

PA 128° 

( 0 ", 0 ") 


9 

F08 

N 2 H+ 

1^0 

93.173 

3.6' 

'x2.7" 

PA 94° 

( 0 ", 0 ") 


14 

Dec 2010-May 2011 

H 

(41) a 

92.034 

3.5' 

'x 2 . 8 " 

PA 94° 

( 0 ", 0 ") 


15 

Dec 2010-May 2011 

CH 3 CN 

53—>4o 

91.971 

3.5' 

'x 2 . 8 " 

PA 94° 

( 0 ", 0 ") 


15 

Dec 2010-May 2011 


52—>42 

91.980 

3.5' 

'x 2 . 8 " 

PA 94° 

( 0 ", 0 ") 


15 

Dec 2010-May 2011 


51^4, 

91.985 

3.5' 

'x 2 . 8 " 

PA 94° 

( 0 ", 0 ") 


15 

Dec 2010-May 2011 


5o^4o 

91.987 

3.5' 

'x 2 . 8 " 

PA 94° 

( 0 ", 0 ") 


15 

Dec 2010-May 2011 

H 

(35)0! 

147.047 

2.3' 

'xl.9" 

PA 75° 

(-1-7", -h2.5"), (-7" 

, -2.5") 

13 

Dec 2010-May 2011 

CS 

3^2 

146.969 

2.3' 

'xl.9" 

PA 75° 

(-1-7", -h2.5"), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 

C-C 3 H 2 

3i.2^22,1 

145.089 

2.3' 

'xl.9" 

PA 75° 

i+l", +2.5”), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 

A-CH 3 OH 

3o,o^2o_o 

145.103 

2.3' 

'xl.9" 

PA 75° 

(+1", +2.5”), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 

E-CH 3 OH 

3o,o^2o,o 

145.094 

2.3' 

'xl.9" 

PA 75° 

i+l”, +2.5”), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 


3-1,0—>2-1,0 

145.097 

2.3' 

'xl.9" 

PA 75° 

(+1”, +2.5”), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 

H 2 CO 

2o,2^1o,i 

145.603 

2.3' 

'xl.9" 

PA 75° 

i+l”, +2.5”), i-T 

, -2.5") 

12 

Dec 2010-May 2011 

HC 3 N 

16^15 

145.561 

2.3' 

'xl.9" 

PA 75° 

(+1”, +2.5”), (-7" 

, -2.5") 

12 

Dec 2010-May 2011 


‘ WOl: Weiss et al. (2001b); GOl: Garci'a-Burillo et al. (2001); F08: Fuente et al. (2008) 


interaction between M 82 and its neighbour M 81, and subse¬ 
quent bar-driven evolution. Each episode lasted only for a few 
million years suggesting strong negative feedback effects from 
the starburst activity. 

X-rays and optical observations have shown the existence of 
a biconical outflow of hot gas coming out of the plane in the 
nucleus of M 82 (Bregman et al. 1995, Shopbell et al. 1998). 
It is widely accepted that the driving mechanism of the outflow 
phenomenon in starbursts is linked to the creation of expand¬ 
ing shells of hot gas by supernovae. These hot bubbles blow out 
into the halo, entraining surrounding cold gas and dust at several 
hundreds of kilometers per second. Interferometric observations 
showed that, contrary to most species, the SiO emission is not 
tracing the molecular gas in the galaxy disk. The detection of a 
-500 pc molecular gas chimney and a super-shell in SiO indi¬ 
cates the occurrence of large-scale shocks in the disk-halo inter¬ 
face (Garcla-Burillo et al. 2001). Out-of-plane molecular emis¬ 
sion has also been detected in CO, HCN and HCO* (Weiss et al. 
2001b, Salas et al. 2014). Recent observations of the HCN 1—>0 
and HCO* 1—>0 lines by Salas et al. (2014) showed that >2% of 
the total dense molecular gas is in the outflow formed by the gas 
expelled by the central starburst. 

M 82 is one of the finest examples of how chemistry can 
help to the full understanding of the interstellar medium (ISM) 
of an external galaxy. Our early HCO interferometric map us¬ 
ing the PdBI showed that the M 82 nucleus is a giant photon- 
dominated region (PDR) of ~650 pc size. Furthermore the com¬ 
parison between the HCO and H^^CO* images suggested that 
the PDR chemistry is propagating across the M 82 nucleus. 
Our subsequent 30m and PdBI chemical studies provided fur¬ 
ther support for the existence of giant PDRs in the galaxy disk 
(Garcla-Burillo et al. 2002, Fuente et al. 2005, 2006, 2008). In 
particular it is remarkable the detection of the reactive ion CO^ 
using the 30m telescope (Fuente et al. 2006), the first time in an 
external galaxy, and the interferometric image of the HOC* 1 —>0 


line (Fuente et al. 2008), the first one ever observed. Chemical 
modeling using the Meudon PDR code showed that most of the 
observations (CO+, HOC*, HCO*, CN, HCN, H 3 O*) are well 
explained assuming that about 87% of the mass of the molecular 
gas is forming small (~ 0.02 pc) clouds with a total thickness of 
~5 mag. In these small clouds all the molecular gas is exposed 
to an intense UV field and the entire cloud is a PDR. Further star 
formation is not expected in this PDR component. A small mass 
fraction (~13%) must, however, be located in shielded regions 
to account for the measured [CN]/[HCN] abundance ratio. We 
estimated that these clouds could have large column densities, 
with a total visual extinction of >50 mag. The existence of this 
shielded component is also necessary to explain the detection of 
some species like NH 3 and CH 3 OH (Weiss et al. 2001a, Martin 
et al. 2006, Aladro et al. 2011). 

However, this chemical study was based on single-dish 30m 
observations towards one position (position E in the nomencla¬ 
ture of Fuente et al. 2008) and cannot be extrapolated to the en¬ 
tire galaxy. In addition, all the observed molecules were PDR 
tracers which are not adequate to trace the shielded gas. In this 
paper, we present high spatial resolution (~100 pc) observations 
of a set of molecular lines carried out with the Plateau de Bure 
Interferometer (PdBI). This set of molecules includes PDR trac¬ 
ers like CN and C-C 3 H 2 but also well-known tracers of cold and 
dense gas like N 2 H*, CS and C^^O. We combine the information 
from all these molecules to gain insights into the spatial distri¬ 
bution, physical conditions and chemistry of the molecular gas 
reservoir in the M 82 starburst. 


2. Observations and data reduction 

The interferometric observations of the CN 1 —>0 (113.491 GHz), 
N 2 H* 1^0 (93.173 GHz), CS 3^2 (146.969 GHz), C-C 3 H 2 
3i, 2 ^ 22,1 (145.089 GHz), H 2 CO 2 o, 2 ^ 1 o.i (145.603 GHz) and 
HC 3 N i6—>15 (145.601 GHz) lines were carried out with 
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Fig. 1. Interferometric spectra towards El [offset (+14", 
+5")] in temperature units (K). The velocity has been set 
to visr-'i3Q km s ', based on the data. Vertical lines 

indicate the following frequencies; 91987 (CH 3 CN Sq—> 4o), 
92034 (H(41 )q'), 145561 (HC 3 N 16^15), 145603 (H 2 CO 
2oi^loi), 145089 (C-C 3 H 2 3i 2 ^ 221 ), 145094 (E- 

CH 3 OH 3oo^2oo), 145097 (E-CH 3 OH 3-io^2_io), 

145103 (A-CH 3 OH 3o,o^2o,o), 146969 (CS 3^2), and 
147047 (H(35)a) MHz. 


the IRAM Plateau de Bure Interferometer (PdBI) as part of 
two different projects. The CN 1—>0 line was imaged in May 
and December, 2005, with the antennas arranged in the C 
and D configurations providing almost a circular beam of 
2.46"x2.27" PA 170°. The primary beam of the PdBI at this 
frequency is ~44", enough to cover the whole galactic plane 
at 3mm and only one field was observed. During the observa¬ 
tions we adjusted the spectral correlator to give a contiguous 
bandwidth of 1 GHz with a frequency resolution of 2.5 MHz. 
The continuum maps were generated using the channels free 
of line emission. Then we subtracted the continuum emission 
to produce the spectroscopic maps. The maps have been cor¬ 


rected for the primary beam attenuation (primary beam=44" at 
113.5 GHz). 

In the second project we imaged the N 2 H* 1—>0 
(93.173 MHz), CS 3^2 (146.969 MHz), C-C 3 H 2 3i, 2 ^ 22,1 
(145.089 MHz), H 2 CO 2 o, 2 ^ 1 o,i (145.603 MHz) and HC 3 N 
16—>15 (145.601 MHz) lines. These observations were carried 
out on December, 2010 and May, 2011, with the antennas ar¬ 
ranged in the C configuration. During the 3mm observations, the 
narrow band correlators were placed to cover the N 2 H* 1 —>0 line 
with a frequency resolution of 20 kHz. We used the wide-hand 
correlator WideX to cover a bandwidth of 3.6 GHz in dual polar¬ 
ization with a channel spacing of about 1.95 MHz. This allowed 
us to observe the H(41)a recombination line at 92.034 GHz and 
the CH 3 CN k=0,1,2,3 at -91.987 GHz simultanously 

with the N 2 H* 1—>0 line (see Pig, 1). Only one field was ob¬ 
served and the maps have been corrected for primary beam atten¬ 
uation (primary beam=54" at 92.5 GHz). The synthesized beam 
was almost circular with HPBW~ 3.56"x2.74" PA 94°. Because 
of the slightly better sensitivity, we used the WideX cube for our 
analysis. 

At 2mm, we observed two fields at (+7", +2.5") and (- 
7", -2.5"). The CS 3—>2 line was covered with the narrow 
band correlators providing a frequency resolution of 160 kHz. 
The wide-hand correlator WideX was used to cover a band¬ 
width of 3.6 GHz in dual polarization with a channel spacing 
of about 1.95 MHz. This allowed us to observe the H(35 )q; 
and the C-C 3 H 2 3 1 ^ 2 —^ 22 , 1 , CH 3 OH 3—>2, H 2 CO 2 q 2 —>lo,i 
and HC 3 N 16—>15 lines in the same setting. In Pig. [1] we 
present the observed spectrum towards the position (+14",+5"). 
Unfortunately, the C-C 3 H 2 line is blended with several CH 3 OH 
3^2 lines and cannot be separated (see Pig. [TJ. The H 2 CO 
2 o, 2 ^ 1 o,i (145.603 MHz) and HC 3 N 16^15 (145.601 MHz) 
lines are partially blended and cannot be imaged separately, 
but can be separated in the outer part of the galaxy where the 
linewidths are narrower. We generated continuum maps using 
the channels free of line emission to subtract the continuum 
from the spectral images. The HPBW of the synthesized beam is 
~2.3"xl.88"PA75°. 

Data analysis and image processing have 
been done using the GIPDAS package software 
(http://www.iram.fr/IRAMPR/GIPDAS/) and JPARSEC 
(http;//conga.oan.es/~alonso/doku.php?id=jparsec). 

All the maps are centered at, RA = 09;55;51.9, Dec = 
+69:40:47.10 (J2000), that is the 2.2 /nm peak as determined 
by Joy, Pester & Harvey (1987). This is also the central posi¬ 
tion used in previous interferometric studies (Weiss et al. 2001b, 
Garcla-Burillo et al. 2001, 2002, Puente et al. 2008). Pester et 
al. (1990) established the 2 yum peak at RA = 09:55:52.4, Dec = 
+69:40:46.00 (J2000), i.e. +2.7" E 1.3" S from our center posi¬ 
tion. The M 82 galactic plane is viewed almost edge-on (i~80°, 
Greve 2011). In order to guide the discussion, we have defined 
four positions El, E2, W2, and W1 [(+14",+5"), (+5",+2"), 
(-5",-2") and (-14",-5"), respectively] that are aligned with 
the center at an angle of 68 °. Within the errors, this is the posi¬ 
tion angle of the galactic plane. The intense continuum source 
SN 41.9+58 is located (-1.25",-1.1") from our position W2 
(Kronberg et al. 1981, Weiss et al. 1999). Position El coincides 
with position E in Puente et al. (2008). The stellar bar is inclined 
-4° relative to the galactic plane, with the western side above 
the adopted major axis of the disk. This explains the tilted ap¬ 
pearance in the molecular emission (see Pigs.and |3]). 

In addition to these new observations, we use previous inter¬ 
ferometric data published by Garcia-Burillo et al. (2001), Puente 
et al. (2008) and Weiss et al. (2001b). In Table 1, we give a brief 


3 





















D. Ginard et al.: The chemical footprint of the star formation feedback in M 82 on scales of '-100 pc** 



Fig. 2. Map of the continuum emission at 3mm (Fuente et al. 2008) and line integrated intensity maps of the 12.8 j-im line of [Nell] 
(Achtermann & Lacy 1995), H(41)a, CO 2—>1, CN 1—>0 and N 2 H* 1—»0 lines. The beam is drawn in the left-down corner of the 
panels and positions El, E2, W2 and W1 are indicated by filled polygons. The position of SN 41.9-1-58 is marked with a white cross. 
Eirst contour is at 3xcr level. Contours levels are: 1.7 mJy/beam, 5 mJy/beam to 35 mJy/beam in steps of 5 mJy/beam (cont-3mm); 
0.13 Jy/beam to 1.04 Jy/beam in steps of 0.13 Jy/beam ([Nell]); 4 to 28 in steps of 4 K x km s“' (H(41) a); 200, 400, 600, 800, 
1200, 1600, 2000, 2400 K x km s^' (CO); 15.8 to 252.8 in steps of 15.8 K x km s^' (CN); 2.5, 5.0, 7.5, 10.0, 12.5 K x km s^' 
(N 2 H*). Colour scale is adjusted to the minimum and maximumm values in each panel. 


^ 10 

a 

o 

o 

m 0 

o 

s 

:§-io 


20 0 -20 

A<x (arcseconds) 



Fig.3. The same as Fig. 2 for the H'^CO* 1^0, HCO 1^0 F=2^1, HOC* 1^0, C 3 H 2 3i, 2 ^ 22 , 1 , C'^O 1^0 and CS 3^2 lines. 
Contour levels are: 0.6 to 4.6 in steps of 0.6 K x km s ' (H'^CO*); 0.7 to 1.9 in steps of 0.2 K x km s ' (HCO); 0.28, 0.70, 
1.12,1.54, 1.96, 2.38 K x km s^' (HOC*); 6.4 to 44.8 by steps of 6.4 K x km s^' (C 3 H 2 ); 2, 8 , 14, 20, 26, 32 K x km s^' (C'^O); 
30 to 150 in steps of 30 K x km s^' (CS). 
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El E2 W2 W1 (mJvAeani) 



-20 -10 0 10 20 
Distance (arcseo) 


Fig. 4. Top: In colour scale, the position-velocity (p-v) diagram 
along the galactic plane (crossing El, E2, W2 and Wl) of 
H(41 )q'. Solid contours are: 2.1 mJy/beam to 16.6 mJy/beam 
by 4.2 mJy/beam. In dashed lines, the same contours but cen¬ 
tered at the frequency of CH 3 CN 5o^4o line (91.987 GHz). 
Bottom: In colour scale, the P-V diagram of the H 2 CO 2o,i^lo,i 
line. Eor comparison, we have drawn in solid contours the lev¬ 
els of the CN 1—>0 line: 10.0 mJy/beam to 86.6 mJy/beam by 
20.0 mJy/beam. In dashed lines, the same contours at the fre¬ 
quency of the HC 3 N 16^15 line. Note that the H 2 CO 2o,i^lo.i) 
and HC 3 N 16—>15 lines are heavily blended in the inner region. 
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description of the observations and their more relevant observa¬ 
tional parameters. 

3. Results 

Fig- [U shows the wide band spectra towards El. The velocity in¬ 
tegrated intensity maps for all the lines detected are shown in 
Eigs.|2]and|3] The maps of the H 2 CO 2 o, 2 ^ 1 o,i (145.603 MHz) 
and HC 3 N 16—>15 (145.601 MHz) are not shown because the 
lines cannot be deblended in the inner region. In Eig. |4l we show 
the position-velocity (p-v) diagrams along the galactic plane 
(straight line across El, E2, W2 and Wl) of the H 2 CO 2 o, 2 —>lo,i, 
HC 3 N 16^15, CH 3 CN and H(41)a lines. The CH 3 CN 

5i—>4i is only detected towards the outer part of the galaxy and, 
at 3cr level, around W2. We subtracted the CH 3 CN 5^—>4^^ line 
from the data cube to create the H(41)a map shown in Eig.|2] In 
Fig. |3] we show the C-C 3 H 2 map although we are aware that the 
line is contaminated with CH 3 OH emission. We will discuss this 
problem in Sect. 4. 

3.1. Galactic plane 

The emissions of all the species studied in this paper are concen¬ 
trated in the galactic plane. However, there are some differences 
in their spatial distributions that we below discuss. 

In Eig.|2]we show the integrated intensity maps of the lines 
that present intense emission towards the inner positions, E2 and 
W2. E2 and W2 are related with the ring of ionized gas that 
is associated with the most recent starburst (~5 Myr) (Eorster- 
Schreiber et al. 2011), This ring is asymmetric with the west¬ 
ern part more intense than the eastern one. The continuum emis¬ 
sion at 3mm and the ionized line [Nell] are maxima towards W2 
(see Euente et al. 2008 and Eig.|2]i. As expected, the emission of 
H(41 )q' follows the same spatial distribution. The emissions of 


Distance (arcaec) 

Fig.5. P-v diagrams along a cut crossing El, E2, W2 and Wl. 
Contours correspond to the line intensity of H(41)a (2.1 to 16.6 
in steps of 2.1 mJy/beam). Colour scale shows the following 
line intensity ratios (in T^, units): A) C'^O 1—>0/H(41)a;; B) CN 
1^0/N2H+ 1^0; C) CN 1^0/[CS] 3^2; D) CS 3^2/N2H+ 
1^0. Horizontal lines correspond to velocities: 100, 270 and 
326 km s ', to make easier the comparison with Eig. 7. Line ra¬ 
tios have been calculated with the original images degraded to 
common angular resolution of 3.8" and a velocity resolution of 
12 km s"'. The C^^O 1—>0/H(41)a line ratio can be considered 
as an observational measurement of the gas illumination by UV 
photons. Note that the CN/N 2 H* and CN/CS line ratios increase 
in the highest UV irradiated regions, i.e., minimum values of the 
C'^O 1^0/H(41)aline ratio. 


the molecular lines are, however, more intense in the outer part 
of the galaxy close to positions El and Wl, with the only ex¬ 
ception of CN that presents its emission peak towards the inner 
ring following the spatial distribution of the ionized lines. CN is 
known to present higher abundances in regions with enhanced 
UV fields (Euente et al. 1993, 1995; Bachiller et al. 1997; Boger 
& Sternberg 2005). Its peculiar spatial distribution suggests that 
the UV radiation from the recently formed stars has a strong im¬ 
pact on the molecular gas chemistry in the inner ring. 

The emission of most molecular lines, H'^CO* 1—>0, HOC* 
1^0, C'^O 1^0, HCO 1^0, C 3 H 2 3 1 , 2 ^ 22 , 1 , and CS 3^2 is, 
however, more intense towards the outer positions, El and Wl 
(see Eig. |3]l. Although all these molecules are brighter in the 
outer part of the galaxy, the exact positions of their emission 
peaks differ, especially in the western part. The C 3 H 2 3i 2 —^3.2 1 , 
CS 3^2, H'^CO* 1^0, and HOC* 1^0 peaks are located 2"- 
4" closer to the dynamical center of the galaxy than HCO 1 ^0. 
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1 

l-0/H(41)a 


Fig. 6. Correlations between the C'^O 1^0/H(41)a, CN 
1—> 0 /N 2 H* 1—>0, and CS 3—> 2 /N 2 H* 1—>0 brighness tempera¬ 
ture line ratios along the galactic plane shown in Fig. 5. The plots 
are in logarithmic scale. In order to get a more reliable result 
avoiding the line borders, we only considered the points within 
the following range of values: 0.5 and 5 for C'^O 1—»0/H(41)a 
line; 5 and 15 for CN 1 ^ 0 /N 2 H+ 1^0; 1 and 6 for CN 1^0/CS 
3—>2 between 1 and 6 . The results of the least square fitting and 
the values of the Pearson correlation coefficients are shown in 
the top and medium panels. 


The intense peak of the C-C 3 H 2 line towards El is very likely 
due to the CH 3 OH contamination (see Sect. 4). 

In Fig. |5] we show the p-v diagrams of the CN 1 ^ 0 /N 2 H* 
1^0, CN 1^0/CS 3^2 and the CS 3 ^ 2 /N 2 H+ 1^0 line ratios 
along the plane defined by positions El, E2, W1 and W2. To 
perform these diagrams, we have degraded the spatial and spec¬ 
tral resolutions of the three line intensity cubes to common val¬ 
ues of 3.8" and 12 km s“*, respectively. We adopted 3xcr as the 
threshold for line detection. There is a clear gradient in the CN 
1 —» 0 /N 2 H* 1 —»0 line ratio along the galaxy plane, being a factor 
of ~3 higher in the inner region than in El. For comparison, we 
show the p-v diagram of the C'^O 1—»0 /H(41)a line ratio that 
can be considered as an observational tracer of the gas illumina¬ 
tion. There are two regions where the CN 1 ^ 0 /N 2 H* 1 ^0 takes 
the peak value of ~15: (i) W2, it is the peak in the H(41 )q; line 
emission and the region with the highest concentration of HII re¬ 
gions; and (ii) E2, the emission of the H(41 )q; line is lower than 
in W2 but the molecular gas column density is also lower, espe¬ 
cially in the ~326 km s“* component that is not detected in N 2 H* 
(see Sect. 3.2); as a consequence, the molecular gas is very likely 
immersed in a high UV field (see Fig. |5l A). The CN 1—>0/CS 
3—»2 line ratio has a ratio ~l-2 along the cut and increases to 


~5 only in the two enhanced UV field regions described above 
(see Fig. |5] C). The variations of these two molecular line ra¬ 
tios are not hazardous but seems related with the distribution 
of the ionized gas suggesting that UV photons are the driving 
agent of the molecular chemistry. In Fig.| 6 ]we show the log-log 
correlation diagrams between the CN 1—» 0 /N 2 H* 1^0 and CN 
1^0/CS 3^2 line ratios, and C'^O 1—»0 /H(41)a. As expected 
there is a good (anti-)correlation between these ratios at low ex¬ 
tinctions (low values of C'*0 1 —»0 /H(41 )q;). At high extinc¬ 
tions, the dispersion is larger since we can have PDRs (diffuse 
clouds, the external layers of giant molecular clouds, dense gas 
around Herbig Ae/Be stars) that do not emit in H(41 )q' but also 
present enhanced CN abundances (see e.g. Fuente et al. 1993, 
1995, Fiszt et al. 2001). On the contrary, although with a large 
dispersion, the CS 3—> 2 /N 2 H* 1—»0 line ratio is ~5 towards all 
positions. 

3.2. Vertical distribution of the molecular emission 

M 82 is associated with a kpc-scale outflow with the gas com¬ 
ing towards us in the northern part. Optical and near infrared 
lines look asymmetric with the most intense half towards the 
north where the radio continuum chimneys are located (Wills 
et al. 1999). This asymmetry is also detected in the emission 
of H(41 )q' which proves that it is not a consequence of the 
dust extinction but to the spatial distribution of the massive star 
forming regions located in the northern half of the disk (see 
Fig.0. On the contrary, the molecular emission extends to ~7.5" 
south from the galactic plane following the molecular super¬ 
shell around the supernova remnant SN 41.9-1-58. Regarding the 
molecular emission, there are important differences between the 
vertical spatial distribution of the different species. 

One interesting case is the comparison of the HCO* and 
HOC* emissions. Fuente et al. (2008) found a clear north-south 
gradient in the [HCO*]/[HOC*] ratio being this ratio maxi¬ 
mum towards the southern part of the supershell associated to 
SN 41.9-1-58. The maximum in the HOC* emission is shifted 
~2.5" north relative to H^^CO*, i.e. closer to the ionized layer 
traced by H(41 )a. Fuente et al. (2008) interpreted this gradient as 
a change in the global ionization degree of the molecular clouds. 

We used the CN, N 2 H*, CS and H(41)a' images to make a 
detailed study of the kinematics and spatial distribution of the 
molecular gas in the direction perpendicular to the galaxy plane. 
In Fig. |7] we show the p-v diagrams in the vertical direction 
across E2, the dynamic center of the galaxy, and W2. In these 
diagrams, the velocity axis has been re-binned to channels of 
12 km s“'. 

Towards E2, several velocity components are detected in the 
H(41 )q' and molecular lines. All these components are within 
the range of velocities of the x2 bar orbits, +120 km s ' (Greve 
2011). The component at ~260 km s * is the most intense 
in molecular emission and is detected in CN, N 2 H* and CS. 
However, the component at ~326 km s * is well detected in the 
CN 1—»0 line, only tentatively detected in CS 3—>2 (3xcr level) 
and remains undetected in the N 2 H* 1 —>0 line. This component 
is associated with a secondary peak in the H(41 )q' emission. In 
the right panel, we show the CN 1—»0 /N 2 H* 1 —»0 line inten¬ 
sity ratio. It varies between ~10 and ~15 along this cut with the 
peaks following the peaks of the H(41 )q' emission. 

Several velocity molecular components are also detected to¬ 
wards W2. In this case, we can also observe a layered structure 
in the vertical direction with the H(41)a and CN emission being 
extended towards north while N 2 H* emission is more extended 
towards the south. Again, there are important chemical differ- 
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Fig. 7. P-v diagrams along cuts perpendicular to the galactic plane at E2 (top), the (0,0) position (middle) and W2 (bottom). Colour 
scale and black contours correspond to the H(41 )q' emission. White contours are intensities of the CN 1—>0, CS 3^2 andNaH^ 1—>0 
lines. First contours (3xcr) and steps are: 10 mJy/beam for CN and CS; 2 mJy/beam for N 2 H^ and H(41)a. The CN 1 ^ 0 /N 2 H^ 1^0 
line intensity (in units) ratio along these cuts is shown in the right column. The images have been degraded to a angular resolution 
of 3.8" and a velocity resolution of 12 km s“^ to perform this line ratio. Vertical lines indicate of the position of the velocity 
components discussed in the text. 


ences between the ~ 100 km s“' (north) and~160kms“^ (south) 
components. The north component is more intense in CN and 
CS, while the south component is brighter in the N 2 H^ emission. 
The CN 1—>0 /N 2 H^ 1 —>0 line intensity ratio shows an excellent 
correlation with the H(41)a emission, with a contrast of a fac¬ 
tor of ~5 between the north and south components. The lowest 
values, ~2-3, are found towards the south, and are even lower 
than those observed in the outer part (E1,W1) of the galaxy (see 
Fig.Ell. 


For comparison, we also show the p-v diagram across the 
dynamic center of the galaxy. The H(41)a emission is detected 
shifted towards the north compared with the molecular emission. 
Similarly to the vertical cut across E2, the CN 1—>0 /N 2 H^ 1 —>0 
line intensity ratio vary between ~10 and ~15 with the peaks 
towards the peaks of the H(41)a emission, suggesting that it is a 
general trend all over the galactic nucleus. 


4. Column densities towards selected positions 

4.1. CN, N 2 H+, C-C 3 H 2 and CS 

The high angular and spectral resolutions provided by our new 
data: CN 1^0, N 2 H+ 1^0, H(41)a, CS 3^2 and C 3 H 2 3^2, 
allow us to carry out a detailed chemical study of the interstellar 
medium in M 82. The spectral cubes of all these lines have been 
degraded to a common angular resolution of 3.8" (ss70 pc) and 
a channel width of 12.8 km s ', that are typical values for the 
size and linewidth of a giant molecular cloud (see e.g. Murray 
2011). From these data, we have selected 12 molecular position- 
velocity peaks (see Fig. [ 8 ]l to carry out our chemical modeling. 
The selection criteria are: the knots must be (i) intense and (ii) 
well identified in velocity and position. In addition, the whole 
set must represent the variety of physical and chemical condi¬ 
tions in the galaxy. The intense and compact emission towards 
these positions minimizes the uncertainties due to spatial filter¬ 
ing effects. The selected knots are listed Table 2 and indicated in 
Fig. | 8 ] and Figs. lA.ll to lA.3l 
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Fig. 8 . Spectral maps of the N 2 H* 1—>0 line convolved to an angular resolution of 3.8". To increase the S/N ratio, the velocity 
resolution has been degraded to 12.8 km s"'. The number in the right-upper corner indicates the central channel velocity. First 
contour and step are 0.028 K (~3xcr). Crosses indicate the dynamical center of the galaxy and the position of SN 41.9H-58. 


GPl to GP 6 are compact molecular peaks along the galactic 
plane that follow the velocity pattern of the xl and x2 bar or¬ 
bits. These knots have been detected in all molecular lines (see 
Fig. | 8 ]and Figs. IA.lllA.3l l. We have added 5 positions, SSI to 
SS5, placed in the supershell associated to SN 41.9-1-58. Note 
that SSI and SS2 are very close to GP4 in the plane of sky 
and are only distinguishable because of their different veloci¬ 
ties. Examining the CN spectral maps, we identified one well 
defined CN emission peak with weak counterpart in the other 
species. We will refer to this position as CN 1 and corresponds to 
the ~326 km s“* component towards E2 (see Eigs.|7]and lA.ll i. 
Eig. lA.3l shows the H(41 )q' PdBI spectral maps. This line probes 
the spatial distribution of the ionized gas that is closely related 
to the local UV field. GP4 is the most intense position in this ra¬ 
dio recombination line. Secondary peaks are spatially coincident 
with GPl, CNl and GPS showing the presence of ionized gas at 
these positions. The absence of H(41)a; emission towards SS2 
and SS3 proves the existence of a lower UV field in the sourth- 
ern part of the supershell. Therefore, our selected positions are 
representative of the different physical conditions, in particular 
the mean UV field, within this galaxy. 


We have estimated the CN, N 2 H*, CS and C 3 H 2 molecu¬ 
lar column densities using the LVG code MADEX (Cernicharo 
2012). Our observations do not allow a multitransitional study. 
Instead, we need to assume a uniform layer with constant physi¬ 
cal conditions. On basis of CO and its isotopologues, Weiss et al. 
( 2001 b) derived densities of n(Fl 2 )~ 5 xlO"^ cm“^ and gas kinetic 
temperatures ranging from aiSO K towards the El and W1 posi¬ 
tions to «il50 K towards E2 and W2. Later, Euente et al. (2005) 
determined densities of 1 x 10 ^ cm“^ from single-dish data of 
the CN 1^0 and 2—»1 lines. Higher densities, ~5xl0^ cm“^, 
were estimated by Euente et al. (2008) using the high excitation 
J=3—>2 and 4^3 lines of HCO*. Bayet et al. (2009) proposed 
two components to fit the CS lines with densities of ~10^ cm^^ 
and 6 x 10 ^ cm~^ for the diffuse and dense components, respec¬ 
tively. These estimates proved that, as expected, different density 
components coexist within our beam. Since we are dealing with 
molecules with large dipole moments, we considered that the 
emission is dominated by dense gas with n(H 2 )>lxlO^ cm“^. 
To estimate the molecular column densities and the uncertain¬ 
ties in our calculations we run a grid of LVG models with n(H 2 ) 
varying from 1x10^-5x10^ cm“^ and Tic-50, 150 K and adopted 
the average value. The hyperfine structures of CN and N 2 H* 
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have not been considered in our calculations because of the large 
linewidths. The estimated column densities are shown in Table 2. 
The errors correspond to the minimum and maximum values ob¬ 
tained in our grid, i.e., these errors do not account for the ob¬ 
servational errors, which are lower, but for the uncertainty in the 
physical conditions. In general, the uncertainties in the obtained 
column densities are within a factor of 2 . 

In our calculations we assumed that all the emission at 
145.089 GHz comes from the C-C 3 H 2 carrier that is not true for 
some positions. Methanol was firstly detected by Martin et al. 
(2006) in M 82. Based on the kinematical and spatial distribu¬ 
tion, they concluded that the emission is mainly coming from 
the intense east and west knots at the outer part of the galaxy. 
Later, Aladro et al. (2011) performed a single-dish multitransi- 
tional study of CH 3 OH anc C-C 3 H 2 towards a position close to 
El and derived their column densities and rotational tempera¬ 
tures. We assumed these values to predict the expected CH 3 OH 
and C-C 3 H 2 line intensities and obtained that ~50% of the emis¬ 
sion must belong to each carrier towards El. According to these 
results the C-C 3 H 2 column density is overestimated by a factor of 
~2 in El. Taking into account the similarties between the molec¬ 
ular chemistry in El and Wl, we consider that, very likely, we 
have a similar contamination, i.e. ~5Q%, towards Wl. 

In Table 2, we also show the [CN]/[N 2 H+], [CS]/[N 2 H+] and 
[CN]/[CS] column density ratios towards the selected molecular 
knots. These values have been estimated using the same proce¬ 
dure as for the column densities. We run a grid of LVG models 
and adopted the average value of the fitted column density ra¬ 
tios. When molecules with similar excitation conditions are se¬ 
lected, this method minimizes the uncertainty due to the assumed 
physical conditions, being the uncertainties in the column den¬ 
sity ratios lower than those of the individual column densities 
(see Table 2). 

Significant variations, i.e. higher than the uncertainties, are 
found in the [CN]/[N 2 H*] ratio across the galaxy. This ratio 
takes values «;30 in GPl and GP 6 located in the outer xl or¬ 
bits and increases to «i80 towards GP4. Values >50 are found 
towards GP5 and SSI, both close to W2. The minimum value, 
[CN]/[N2 H*]ss 19, is found in the southern part of the supershell. 
In Fig.|9]we compare the obtained column density ratios with the 
integrated intensity maps of the N 2 H* 1—>0, H(41)a' lines and 
SiO 2—>1 lines. It is remarkable that there is a good correlation 
between the values of the [CN]/[N 2 H*] ratio and the emission 
of the H(41)a line, supporting the interpretation that the vari¬ 
ations of the [CN]/[N 2 H*] ratio are related to the variations of 
the local UV field. We do not find any correlation between the 
[CN]/[N 2 H*] ratio and the SiO emission. 

The [CN]/[CS] ratio behaves similarly to [CN]/[N 2 H+]. It 
takes values *3.5 in GPl and GP 6 located in the outer xl orbits 
and increases to >5 towards GP4, GP5, and SSI. In the case of 
CNl, we estimate a higher value of, [CN]/[CS]=sl0. There is a 
good correlation between the [CN]/[CS] ratio and the H(41)a 
emission. 

Contrary to the [CN]/[N 2 H+] and [CN]/[CS] ratios, the 
[CS]/[N 2 H*] ratio remains quite constant across the galactic 
plane with values around ~12. There is no evidence of varia¬ 
tion of this ratio with the local UV field. When comparing the 
galactic plane positions with those of the supershell, we realized 
that the ratio tends to be a factor of 2 lower in the southern part 
of the supershell, although always within the uncertainty of our 
calculations. 

We used the C'^O map from Weiss et al. (2001b) to deter¬ 
mine the total molecular hydrogen column densities and abso¬ 
lute fractional abundances. The C'^O column densities vary by 


only a factor of ~2 among the different knots with the largest 
values in GPl and GP 6 . Molecular abundances have been de¬ 
rived from the N(X)/N(C'®0) column density ratio assuming 
[C1®0]/['2 cO]= 0.005 and ['2CO]=4xlO-^ which are the aver¬ 
age values derived by Weiss et al. (2001b) on basis of their mul- 
titrasitional study of CO, '^CO and C'^O. Note that these val¬ 
ues correspond to a ratio of 200 which is different from 

the Milky Way standard value of ~500-600. We are assuming 
a constant C'^O abundance all over the galaxy which is an ap¬ 
proximation since the C^^O abundance is dependent on the en¬ 
vironment (UV flux, density, temperature). However, we are in¬ 
terested in the dense molecular gas (n(H 2 )>lxlO^ cm“^) and the 
C'^O abundance is expected to remain roughly constant in this 
component. Following this procedure we derived the molecular 
abundances shown in Table 2. The N 2 H* abundance is ~3xl0“^' 
for all the positions, corroborating that N 2 H* is a good tracer of 
the dense molecular gas. An abundance of a few 10“^° is mea¬ 
sured for CS. Small changes, less of a factor of 2, are found in 
the abundances of this species that are within the uncertainties 
of our calculations. Only the CN abundance presents significant 
variations within the galaxy, being larger by a factor ~3 in GP4 
than in the outer posisions, GPl and GP 6 , and the southern su¬ 
pershell positions, SS2 and SS3. 

4.2. Other species: H 2 CO, HC 3 N and CH 3 CN 

In order to have a more comprehensive view of the chemistry in 
M 82 we have compared our new data with previous interfero¬ 
metric images. For this aim, we degraded the angular resolution 
of all our images to ss5.9” (ail 11 pc), which is the lowest from 
the data set shown in Table 1, and used the total velocity inte¬ 
grated line intensity in our calculations. The Gaussian fit param¬ 
eters towards El, E2, W2 and Wl are shown in Table A.l. 

Column densities and column density ratios were calculated 
following the same procedure as in Sect. 4.1. In Table 3, we show 
the values obtained for El, E2,W2 and Wl. The [CS]/[H'3 cO*] 
and [H'^C0^]/[C*^0] are quite uniform across the galaxy. This 
suggests that CS, H'^CO* and C'*0 are good tracers of dense 
gas in a wide range of physical conditions. According to these 
results, El and Wl present the highest total molecular hydro¬ 
gen column densities. This is consistent with the detection of 
complex molecules like NH 3 , HC 3 N, CH 3 CN and CH 3 OH to¬ 
wards them (Weiss et al. 2001a, Martin et al. 2006, Aladro et al. 
2011). The small hydrocarbons C-C 3 H 2 and C 4 H are also good 
tracers of low and high UV PDRs (Fuente et al. 2003, Pety et al. 
2005, Pilleri et al. 2013, Cuadrado et al. 2014). Taking into ac¬ 
count that the abundance of C-C 3 H 2 is very likely overestimated 
by a factor of ~2 in El and Wl (see Sect. 4), we conclude that 
the abundance of C-C 3 H 2 could be enhanced by a factor of ~3 
in the inner x 2 orbits, in agreement with the trend observed in 
CN. Moreover, the derived H 2 CO abundance is similar to that 
derived by Guzman et al. (2013, 2014) in the Horsehead neb¬ 
ula and consistent with the interpretation of the H 2 CO emission 
coming from PDRs. Summarizing, the overall chemical behav¬ 
ior observed in M 82 is well interpreted as the consequence of 
the effect of UV radiation on the molecular chemistry. This pro¬ 
duces an increase in the abundance of PDR tracers (CN, C-C 3 H 2 , 
H 2 CO) towards E2 and W2, following the spatial distribution of 
HII regions. The detection of complex molecules shows that a 
fraction of the gas is protected from the UV radiation in the inte¬ 
rior of large molecular clouds. The amount of gas in this shielded 
component is higher in the outer xl bar orbits. 

We remind that the position W2 is related with the super¬ 
nova explosion SN 41.9H-58 that is launching ionized gas out 
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Fig. 9. [CN]/[N 2 H*], [CS]/[N 2 H*] and [CN]/[CS] abundance ratios for the positions listed in Table 2 in the galacic plane (left) 
and the supershell associated with SN 41.9+58 (right) superposed to the N 2 H* 1—>0 integrated intensity map (grey scale). For 
comparison, on the left column we plot the contours of the integrated intensity emission of the H (41) cr line, and on the right, the 
integrated intensity emission of the SiO 2—»1 line (Garcfa-Burillo et al. 2001). The white cross indicates the position of SN 41.9+58. 


Table 2. Selected regions for the chemical model'. 


Point 

Offset 

Vel 

N(C'“0) 

X(CN) 

XfNjH*) 

Xfo-CsHa)^ 

X(CS) 

[CN]/[N2H*] 

[CSMNjH*] 

[CN]/[CS] 


(") 


(km s“') 

(xlO'® cm“^) 

(xlO-®) 

(xlO-") 

(xl0-'°) 

(xlO-'") 




GPl 

GP2 

+ 12.1 
+ 11.6 

+4.1 
+ 1.0 

327.4 

289.0 

2 2*'-'' 

0 9+u» 

10^8:^ 

^•^-0.3 

3 O*' *' 

o '-14-3.0 

26*i:^ 

^•^-1.4 

3 0*“’ 

-1.3 

33*" 

^^-13 

12+^ 

^ -6 

13)1° 

+(Tq3 

0 0 4-0.5 
'^•'^-0.5 

GP3 

+6.6 

-0.8 

263.4 

1 8+10 

i-n" 


3-0)i'5 

3.8!“ 

49+23 

^^-16 

14)10 

Q 74-O.6 
-0.7 

GP4 

-5.1 

-2.0 

97.0 

0 9+0-5 

3-O-u 

4 3+2-0 

2 6+2-« 

3 5+^-2 

80)^^ 

11*^ 

7.4)[^ 

4 8*°:° 

GP5 

-10.4 

-4.7 

97.0 

1 9+^-^ 

^■^-0.9 


3 2+^-^ 

2 4+2-'7 

2 

^•^-1.0 

49+23 

^^-16 

11)1 

GP6 

-21.1 

-5.6 

84.2 

2 l + 

0.8)“ 



1 8+2-5 

28)12 

7)3 

'1 04-0.6 
■+‘°-0.6 

551 

552 

553 

554 

-1.2 

-2.6 

-9.4 

-13.3 

-2.0 

-7.8 

-7.6 

-2.9 

173.8 

161.0 

148.2 

148.2 

1 0*°'^ 
06*8:^ 

g78:? 

2 O*' ® 
09^8:2 

0 9^^^ 

4 8*^'^ 

27*1:^ 

2771:1 

' 13 

3.1)“ 

1 4+l:s 

22*8:? 

Z,.Z, . Q 

9 q+4.3 
11*1:8 
07?:? 

1 1 

47*23 

19)” 

38*11 

66)i 

4*3 

5*^ 

5+4 

6 2*'-° 
75*1:^ 

427A:t 

40*8:8 

-9-7 

8.3);;° 

SS5 

-2.8 

-2.5 

135.4 

1 

^'^-0.5 

^•^-0.8 

A 9+^-f^ 
^•^-2.0 

^•^-1.3 


10)J 

CNl 

+4.0 

+2.0 

327.4 

-Q.l 

2 0**'* 

_ ^•''-0.7_ 

<4.0* 

<3.0* 


>50 

>4 

9 g+2.0 
^-^-2.0 


'Relative abundaces wrt. H2; ^ Contaminated with methanol (see text); * 3cr limit. 


of the galaxy plane. The detection of SiO by Garcfa-Burillo et 
al. ( 2001 ) proved the existence of shocks and their effect on the 
gas chemistry. H 2 CO, CH 3 OH and CS are considered good trac¬ 
ers of shocks in galactic and extra-galactic environments (see 
e.g. Bachiller et al. 2001, Garcfa-Burillo et al. 2000, Usero et 


al. 2006). However, our data do not present any evidence of CS 
abundance enhancement because of possible shocks. In fact, the 
[CS]/[N 2 H*] ratio is slightly lower in the super-shell. This sug¬ 
gests that only a small fraction of the gas is affected by shocks, 
which is also consistent with the low average SiO abundance 
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Table 3. Fractional abudances and abundance ratios. 



El 

E2 

W2 

W1 

N(C‘’*0) (X 10'") 

10!° 

4 1+^-5 

^•°-2.0 

ll!i 

X(CN) (X lO-®) 

1 4+0.5 

^ -^-0.4 

Q 9 + 1.0 
-0.8 

-1.0 

l-2!^i 

X(N2H*) (X 10“") 

9 4+0.6 
0.6 

■^•^-0.4 

4 4+0.8 

^•^-0.9 

^•^-0.4 

X(CS) (X 10-“*) 

5 

*^•^-1.8 

*7 '7 + 3.6 

2.5 

*7 q+3.6 
' •°-2.8 

4 

^••^-1.4 

X(H‘-’CO+) (X 10“") 

1 q+0.7 
^•^-0.7 


9 4+10 
-0.8 

13 *°° 

X(o-C3H2) (X 10-'°) 

-1.2 

2-8!?l 

3.4!|° 

-1.6 

X(p-H2CO) (X 10-") 

7-4!|* 

'^^•^-4.8 

>5.6!°:^ 

4.6*?-? 

X(HC3N) (X 10-*') 


<13+19* 

<3i!^r 

5 3+^-^ 
-4-2 

X(CH3CN) (X 10-'*) 


< 1 . 0 ** 

<0.8** 

1 

^•^-0.4 

[CN]/[N2H^] 


101!?; 

92!|« 

58!*^ 

[CA]/[CS] 

2-8!!:? 

4 7+1-5 

^•'-1.7 

5.6!°? 

9 Q+0.9 
^•^-0.9 

[CS]/[N2H*] 

20!° 

17 

^ ' -8 

23!° 

2U? 

[//'"CO+l/lW//*] 

0-2°i2 

^‘^-0.2 

0.5!°? 

0.6!°-, 


Notes. Beam (5.9") average column densities resulting from LVG cal¬ 
culations in a grid with Tj.=50-150 K, and n(H2)=lxlO^, 5x 10^. 
Errors correspond to the minimum and maximum values in the grid. 
* Reasonable guess since the HC3N and H2CO are blended at this 
position.** Upper limit have been derived assuming Av=50 km s“*, 
Tt=50 K and n(H2)=5xl05 cm-T 

determined by Garcia-Burillo et al. (2001), X(SiO)~lxlO 
when compared with the SiO abundances, ~10“^-10’®, found in 
the shocks associated with galactic star forming regions (Martm- 
Pintado et al. 1992, Bachiller et al. 2001). In the case of H 2 CO 
anc CH 3 OH we cannot do a detailed analysis because their lines 
are blended with other species. 

5. Chemical model 

We used the Meudon PDR code 1.4.4 (Le Petit et al. 2006) to 
model the chemistry in M 82. In our calculations, we assume 
that the ISM is composed of clouds bathed by an intense UV 
field. We simulate each cloud by a uniform plane-parallel layer 
illuminated from the two sides. This layer has constant density 
(n=n(H)-i- 2 xn(H 2 )= 4 xlO^ cm“^) and the gas kinetic tempera¬ 
ture is calculated by detailed heating and cooling balance. The 
adopted initial elemental abundances are the same as in Fuente 
et al. (2008). We run a grid of models varying the cloud size, UV 
field and cosmic rays ionization rate. 

Our grid includes two values of the Draine field, 6x10^ and 
6x10^. The first value was determined by Fuente et al. (2005, 
2006, 2008) from the CO+, HCO* and HOC* observations. The 
second value is considered to account for the possibility of re¬ 
gions with a lower UV field within the galaxy. To investigate 
the effect of the enhanced cosmic-ray flux on the chemistry, 
we have repeated our model calculations with ^=5x10 and 
5x10“^^ s"'. The first value corresponds to the cosmic rays ion¬ 
ization rate in the molecular clouds of our Galaxy (Indriolo et 
al. 2014). The second is the value estimated by Suchkov et al. 
(1993) to account for the physical conditions of the molecular 
gas in M 82. We are aware that X-rays are intense and could 
also play an important role in the chemistry (Fuente et al. 2008, 
Spaans & Meijerink 2007). Our code does not include the X- 
rays in the chemical calculations. We investigate the effect that 
X-rays could have on the chemistry by increasing the cosmic 
rays ionization rate to 5x10“'"^ s“'. To increase the cosmic rays 
ionization rate reproduces reasonably well the effect of X-rays 
on the molecular chemistry although cannot account for the ef¬ 


fect on the molecular gas heating. Finally, in our simple model, 
the cloud size is given by the total cloud visual extinction. We 
considered five values; 5 mag, 10 mag, 20 mag, 50 mag and 
100 mag. 

In Fig. [To] and Table 4 we show the cloud average molecu¬ 
lar abundances as a function of the different model parameters. 
The large dispersion in the values is not unexpected since all the 
parameters are varied by several orders of magnitude. The cloud 
average abundances of CN, N 2 H^, CS and HCO* increase with 
the cosmic rays ionization rate. This increase is especially im¬ 
portant in the case of N 2 H*. While the abundances of the other 
species increase by a factor of a few, the abundance of N 2 H* 
increases by two orders of magnitude between ^=5x10“'^ and 
5x10“'"^ s"'. Surprisingly, the N 2 H* abundance is quite constant 
in the studied knots. This suggests that the chemical changes we 
observe are not related with variation in the X-rays flux. We can¬ 
not discard, however, that variations in the UV field and cloud 
sizes could balance to keep the N 2 H* abundance constant in 
spite of the variation of the X-rays flux (Strickland & Heckman 
2007). 

The cloud size is the parameter that has the largest impact 
on the results of our model. N 2 H* has a negligible abundance in 
clouds with Ai,<10 mag and abundances varying between 10“'° 
and 10 ° for larger clouds. This corroborates our assumption that 
this species is an excellent tracer of the UV-protected molecu¬ 
lar gas. The [CN]/[N 2 H*] ratio is a good measure of the cloud 
size in the case of one single cloud component. Fuente et al. 
(2008) proposed the existence of two cloud components towards 
El; most of the mass, ~87%, is locked in small clouds of Ay~5 
mag and the rest is forming A,,>50 mag clouds. In this case, the 
[CN]/[N 2 H*] ratio would depend on both, the fraction of each 
cloud component and the assumed size for the large clouds. 

The cloud average CS abundance is quite constant for sizes 
between ~5 and ~20 mag but increases by almost one order of 
magnitude for larger sizes. The large values of the CS abundance 
in the 50 mag clouds are not reliable. The CS abundance is very 
dependent on the assumed elemental sulphur abundance. Our 
gas phase models does not consider the adsorption of molecules 
on the grain surfaces. Within our galaxy, the sulphur abundance 
measured in PDRs is a factor >4 lower than the solar value, sug¬ 
gesting that significant sulphur depletion occurs even in highly 
irradiated environments (see Goicoechea et al. 2006). Moreover, 
our model assumes the solar value for the sulphur elemental 
abundance, 1.8x10“^, that could greatly differ from that in M 82 
(Umeda et al. 2002, Origlia et al. 2004). Therefore, we do not 
consider that the values of the CS abundance predicted by our 
model are reliable. 

The [CN]/[N 2 H*] ratio is a sensitive tracer of the cloud size 
for sizes <50 mag. Other parameters like the incident UV field 
and the cosmic rays ionization rate could also affect the pre¬ 
dicted [CN]/[N 2 H+] ratio. As expected, the [CN]/[N 2 H+] ratio 
increases for higher UV fluxes, especially for small clouds. On 
the contrary, the [CN]/[N 2 H*] ratio decreases with the cosmic 
rays ionization rate because of the increase of the N 2 H* abun¬ 
dance. In the following, we discuss the cloud distribution in 
M 82 on basis of the [CN]/[N 2 H*] ratio. We use our a priori 
knowledge of the galaxy to constrain the chemical models. We 
discuss the results for one- and two-component models. 

In Fig. [TT] we plot the [CN]/[N 2 H*] abundance ratio as a 
function of the cloud size for Draine fields of 6x10° (high-UV) 
and 6x10° (low-UV). We adopt the value derived by Suchkov 
et al. (1993) for the cosmic rays ionization rate because it is 
the most likely. GP4 is the position with the highest value of 
the [CN]/[N 2 H*] ratio and also the emission peak of the H(41)a 
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Fig. 10. Model predictions as a function of the cosmic ray ionization rate, Draine field and cloud size. Model parameters and results 
are shown in Table 4. 


line. It is reasonable to think that high-UV models, in red lines 
in Fig.[TT] are more adequate to account for the molecular abun¬ 
dances at this position. One single cloud component with sizes 
of ~22 mag could explain the observed [CN]/[N 2 H*] ratio of 
~80. Lower [CN]/[N 2 H*] ratios, ~30, are measured towards the 
positions GPl and GP 6 . Although less intense, these knots are 
also associated with peaks in the H(41)a line emission which 
proves the presence of energetic UV photons capable to ionize 
the atomic hydrogen (see Fig. lA.3l) . Assuming the high-UV case, 
the lower [CN]/[N 2 H*] ratio can be explained by the existence 
of a population of larger clouds, ~30 mag (see Fig. fTTI) . The 
existence of larger clouds is also consistent with the lower ki¬ 
netic temperatures and higher densities measured by Weiss et al. 
( 2001 b) and the detection of complex molecules at these posi¬ 
tions. Low values of the [CN]/[N 2 H*] ratio, ~20-30, are also 
observed towards the SS2 and SS3 points placed in the molecu¬ 
lar supershell. Since there is a lack of H(41) a emission towards 
these positions, we used the low-UV field models to interpret 
the [CN]/[N 2 H*] ratio and obtained that the molecular emission 
could come from translucent clouds of ~14 mag. 

The variation of the average cloud size we conclude from 
our one-component model is compatible with the results of 
Fuente et al. (2008) who interpreted the abundances in El 
in a two-component scenario. We can use a two-component 
model to explain our observations. In this case, the CN emission 


would come from the small (with Av~5 mag) and large (with 
Av~50 mag) molecular clouds while the N 2 H* emission would 
come only from large molecular clouds. Fuente et al. (2008) es¬ 
timated that only ~13% of the molecular gas is forming large 
molecular clouds in EL Assuming the same fraction at all posi¬ 
tions, the [CN]/[N 2 H*] ratio in the large cloud component would 
be about 0.13 times the observed one, i.e. ~10 in GP4, ~4 in GPl 
and ~3 in SS3, i.e., we would have clouds of >50 mag in the 
three positions, with the largest ones, ~100 mag, towards GPL 
Summarizing, both, the one and two-components models, come 
to the conclusion that the highest concentration of UV-protected 
molecular gas is found towards the xl orbits. The detection of 
complex molecules is better understood in terms of two popula¬ 
tions of clouds with a small fraction of the gas, ~13%, locked in 
giant molecular clouds where the gas is UV protected and can 
form large molecules. 

Of course, the two models are very simple and have to be un¬ 
derstood as a guide to interpret the molecular chemistry in this 
starburst galaxy. As commented above, our model neglects the 
surface chemistry and X-rays effets that could be important for 
some species (see discussion above). Moreover our gas phase 
model assumes steady-state chemistry. The characteristic time 
for PDR chemistry is t~1I{Kci Gq), where Kd is the photocissoci- 
ation rate in s ' and Go is the incident UV field in Habing units. 
The value of kj depends on the extinction from the illuminat- 
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ing source as Kd~KoX exp(-b A,,). Assuming a typical value of 
s b=1.8 and the incident mean UV interstellar field in 
M 82, Go=10"^ Habing fields, the characteristic time is >10 Myr 
for Ai,>2.5 mag. This means that if the last starburst episode took 
place about ~5 -10 Myr ago, the chemistry might be out of equi¬ 
librium in the clouds interior. Our steady-state scenario is there¬ 
fore only adequate for the cloud surfaces and small clouds. The 
shortcoming of time-dependent models is that depend on the ini¬ 
tial conditions and require of a more sophisticated analysis that 
is beyond the scope of this paper. Subsequent improvement of 
our current interpretation would require to couple the star for¬ 
mation history in this galaxy with a time-dependent PDR code 
including surface chemistry (Bayet et al. 2009, Viti et al. 2014, 
Bisbas et al. 2012). 

For all species, the estimated fractional abundances in M 82 
are lower than those predicted by our gas-phase model. Absolute 
fractional abundances are difficult to compare with models be¬ 
cause of the large uncertainty in the total molecular gas column 
density. Weiss et al. (2001b) derived the total molecular gas col¬ 
umn densities in M 82 using different methods; LVG calcula¬ 
tions, Local Thermodynamic Equilibrium (LTE) solution and the 
standard Xco conversion factor, and their results agree within 
a factor ~3-4. In this paper, we derived the molecular column 
density towards the different points using the C'^O data and 
assuming the CO abundance derived by Weiss et al. (2001b), 
X(CO)=5xlO“^. We would like to note that this value is already 
lower than the value predicted by our model for large clouds 
in which essentially all the carbon is in CO with an abundance 
of ~1.3xl0 Therefore our fractional abundances are accurate 
within a factor of 4-10 and consequently, roughly consistent 
with model predictions. 

6. PDR ans shock chemistry in M 82 

M 82 has been subject of succesive starburst episodes in the last 
~200 Gyr. The most recent (~5 Myr) took place throughout the 
central regions of M 82 and was particularly intense in the x2 
orbits and along the stellar bar. The negative feedback effects 
from this starburst activity is determining the existence and dura¬ 
tion of future starburst episodes. We have carried out a chemical 
study of the molecular gas at scales of ~ 100 pc in order to inves¬ 
tigate the feedback effects and the future star formation activity 
of the galaxy. 

Our main result is that the chemistry of the molecular gas is 
determined by the intense UV radiation produced by the massive 
stars. There is a systematic variation of the studied molecular 
abundance ratios with the intensity of the H(41)a emission that 
corroborates that the whole nucleus (~650 pc), down to scales 
of -100 pc, behaves like a giant PDR. Within this PDR, we have 
chemical variations due to changes in the local UV field and 
in the structure of the molecular clouds. There are three well 
differentiated regions: 

- The inner x2 bar orbits that are associated with the most re¬ 
cent starburst and where the PDR tracers present their peak 
abundances. The detection of N 2 H* in these regions proves 
that molecular clouds of >20 mag are present and there¬ 
fore, the molecular gas reservoir to form new stars is not 
exhausted. 

- The outer xl bar orbits where most the UV-protected molec¬ 
ular gas is located. The detection of H(41)a shows that mas¬ 
sive stars have already been formed in these clouds. 

- The molecular supershell associated with the supernova rem¬ 
nant SN 41.9-1-58. There are no signs of massive star forma- 




Fig. 11. [CN]/[N 2 H*] abundance ratio as a function of the cloud 
size for models with CR=5x 10 '^ s ' and Draine fields of 6x10^ 
(red line) and xlO^ (blue line). Top Cloud sizes in the one- 
component model. Bottom Sizes of the large cloud component 
in the two-component model (see text). The values and errors for 
GPl, GP4 and SS3 are indicated. 


tion activity in the southern part of this super-shell although 
there are intense N 2 H* knots that could be sites for future 
star formation. 

The detection of a -500 pc molecular gas chimney and a 
super-shell in SiO indicates the occurrence of large-scale shocks 
in the disk-halo interface (Garcfa-Burillo et al. 2001). We re¬ 
mind that the supershell is related with the supernova explosion 
SN 41.9+58 that is launching ionized gas out of the galaxy plane. 
However, our data do not present any evidence of possible shock 
chemistry. This suggests that only a small fraction of the dense 
molecular gas is affected by shocks, which is also consistent 
with the low SiO abundance determined by Garcfa-Burillo et al. 
(2001), X(SiO)-lxlO when compared with the SiO abun¬ 
dances, -10“®-10“®, found in the shocks associated with galac¬ 
tic star forming regions (Martin-Pintado et al. 1992, Bachiller et 
al. 2001). Higher spatial resolution observations are required to 
detect the layer of molecular gas around the supershell whose 
chemistry is dominated by shocks. This is also consistent with 
the moderate fraction of dense molecular gas (-2%) that Salas 
et al.(2014) detected in the halo, in contrast with -25% found 
by Walter et al. (2002) on basis of CO observations. Most of the 
expelled gas presents moderate densities and is not detected in 
dense molecular tracers. 
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Table 4. Chemical models' 


Model 

CR (xlO-'O 

Cloud size (mag) 

Draine field 

X(CN) 

XfNjH*) 

X(CS) 

[CN]/[N2H*] 

[CN]/[CS] 

[CS]/[N2H*] 

1 

5000 

100 

6e3 

High-UV 

1.2e-7 8.0e-9 

2.1e-7 

15 

0.56 

27 

2 

500 

100 

6e3 

1.2e-8 

2.8e-9 

2.4e-7 

4.2 

0.05 

85 

3 

5 

100 

6e3 

1.3e-9 

l.le-10 

4.7e-8 

12 

0.03 

439 

4 

5000 

50 

6e3 

l.le-7 

4.6e-9 

l.Oe-7 

23 

1.0 

22 

5 

500 

50 

6e3 

1.4e-8 

1.2e-9 

7.0e-8 

11 

0.2 

58 

6 

5 

50 

6e3 

2.4e-9 

8.1e-ll 

5.2e-8 

29 

0.05 

641 

7 

5000 

20 

6e3 

l.Oe-7 

8.3e-10 

1.2e-8 

130 

9.0 

14 

8 

500 

20 

6e3 

2.2e-8 

2.2e-10 

5.0e-9 

100 

4.4 

22 

9 

5 

20 

6e3 

l.Oe-8 

2.4e-10 

8.0e-8 

2173 

0.07 

333 

10 

5000 

10 

6e3 

5.7e-8 

5.3e-ll 

1.7e-9 

1055 

38 

32 

11 

500 

10 

6e3 

2.0e-8 

2.1e-ll 

1.6e-9 

952 

12 

76 

12 

5 

10 

6e3 

1.2e-8 

l.le-12 

1.9e-9 

1.0e4 

6 

1727 

13 

5000 

5 

6e3 

1.9e-8 

3.0e-13 

l.Oe-9 

6.4e4 

19 

3333 

14 

500 

5 

6e3 

2.0e-8 

6.4e-14 

1.2e-9 

3.1e5 

16 

18750 

15 

5 

5 

6e3 

2.2e-8 

3.6e-14 

1.4e-9 

6.1e5 

16 

38889 

16 

5000 

100 

6e2 

Low-UV 
l.le-7 l.Oe-8 

3.1e-7 

11 

0.38 

30 

17 

500 

100 

6e2 

8.6e-9 

4.0e-9 

3.0e-7 

2.1 

0.03 

74 

18 

5 

100 

6e2 

9.9e-10 

l.Oe-10 

2.8e-8 

9.4 

0.03 

270 

19 

5000 

50 

6e2 

l.le-7 

8.2e-9 

2.4e-7 

13 

0.45 

29 

20 

500 

50 

6e2 

9.4e-9 

3.2e-9 

2.4e-7 

3 

0.04 

75 

21 

5 

50 

6e2 

1.8e-9 

8.4e-ll 

3.2e-8 

122 

0.06 

380 

22 

5000 

20 

6e2 

l.Oe-7 

5.5e-9 

1.7e-7 

19 

0.59 

31 

23 

500 

20 

6e2 

1.3e-8 

1.8e-9 

1.9e-7 

7 

0.07 

105 

24 

5 

20 

6e2 

4.6e-9 

4.6e-ll 

5.0e-8 

98 

0.01 

1086 

25 

5000 

10 

6e2 

8.9e-8 

1.9e-9 

7.2e-8 

47 

0.81 

37 

26 

500 

10 

6e2 

2.0e-8 

3.8e-10 

7.6e-8 

52 

0.26 

200 

27 

5 

10 

6e2 

9.0e-9 

4.3e-12 

2.0e-8 

28 

0.45 

4651 

28 

5000 

5 

6e2 

4.0e-8 

2.8e-ll 

3.2e-9 

1428 

12 

114 

29 

500 

5 

6e2 

1.8e-8 

1.3e-12 

3.0e-9 

1.4e4 

6 

2307 

30 

5 

5 

6e2 

1.7e-8 

4.0e-15 

3.6e-9 

4.3e6 

5 

9E5 


' Notation: 6e3 = 6x10^ 


Summarizing, we present a comprehensive chemical study 
of the molecular gas in the starburst galaxy M 82 using high spa¬ 
tial resolution (60-100 pc) interferometric images. Our chemical 
study shows that the feedback effects are strong in this evolved 
starburst. In fact, the whole nucleus presents a PDR-like chem¬ 
istry which suggests that UV radiation is driven by the subse¬ 
quent evolution of the ISM. Even though, the detection of N 2 H* 
proves the existence of dense molecular gas enclosed in clouds 
of Ai,>20 mag that could form new stars. The main reservoir 
of the dense molecular gas is located in the outer xl orbits that 
might be the preferred site for a new generation of stars proceed¬ 
ing with the inside-outside scenario. 
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Table A.l. Gaussians parameters to the observed spectra' 


Line 


Freq(GHz) 

El (+14" 
I(Kkms-') 

, +5") 

v,jr(km s^‘) 

Av(km s *) 

Tmb(K) 

CN 

1^0 

113.491 

64.15 (0.56) 

316.26 (0.33) 

75.73 (0.77) 

0.79 

H'^CO* 

1^0 

86.754 

4.52 (0.67) 

324.66 (5.02) 

69.39 (12.4) 

0.06 

H 

41a 

92.034 

5.42 (0.19) 

320.48 (1.21) 

70.79 (2.89) 

0.07 

CH3CN 

5t^4j.* 

91.987 

2.43 (0.20) 

328.58 (3.13) 

74.26 (6.62) 

0.03 

N2H+ 

1^0 

93.173 

6.98 (0.17) 

321.42 (0.73) 

60.23 (1.57) 

0.11 

c'*o 

1^0 

109.782 

22.74 (0.62) 

324.40 (0.93) 

69.58 (2.18) 

0.31 

C3H2 

3^2 

145.089 

9.90 (0.24) 

315.23 (0.75) 

63.77 (1.69) 

0.14 

H2CO 

2 o, 2 —>lo,l 

145.603 

8.45 (0.20) 

321.30 (0.67) 

56.90 (1.58) 

0.14 

HC3N 

16^15 

145.561 

4.35 (0.17) 

330.32 (0.90) 

43.12(1.96) 

0.09 

cs 

3^2 

145.760 

31.84 (0.37) 

322.56 (0.32) 

55.63 (0.71) 

0.54 




E2 (+5", 

+ 2 ") 



Line 


Freq(GHz) 

I(Kkms-') 

v,jr(km s^‘) 

Av(km s“*) 

Tmb(K) 

CN 

1^0 

113.491 

56.09 (0.53) 

249.23 (0.35) 

79.56 (0.92) 

0.66 




4.83 (0.30) 

330.76 (0.78) 

26.26 (1.70) 

0.17 

H'^CO* 

1^0 

86.754 

2.19(0.30) 

258.66 (4.28) 

57.14 (7.48) 

0.03 

H 

41a 

92.034 

7.25 (0.24) 

278.04 (1.89) 

121.01 (5.46) 

0.06 

CH3CN 

5t^4t* 

91.987 

1.35 (0.20) 

284.18 (5.79) 

85.36 (13.7) 

0.01 

N2H+ 

1^0 

93.173 

3.74 (0.15) 

255.80 (1.08) 

58.24 (2.89) 

0.06 

c'*o 

1^0 

109.782 

9.22 (0.65) 

259.28 (1.78) 

53.14(4.66) 

0.16 

C3H2 

3^2 

145.089 

4.95 (0.17) 

245.88 (0.88) 

53.58 (2.13) 

0.09 




0.57 (0.12) 

311.92 (2.90) 

26.65 (5.71) 

0.02 

H2CO 

2o,2“>1o,i 

145.603 

4.34 (0.16) 

256.14(0.77) 

45.57 (1.98) 

0.09 

HC3N+H2C02 

145.561 

2.85 (0.17) 

253.33 (1.48) 

53.47 (3.90) 

0.05 

cs 

3^2 

146.969 

16.83 (0.15) 

258.63 (0.20) 

46.79 (0.53) 

0.34 




1.43 (0.11) 

329.47 (0.94) 

25.13 (2.28) 

0.05 




W2 (-5" 

, - 2 ") 



Line 


Freq(GHz) 

I(Kkms-') 

v,jr(kms^') 

Av(km s“') 

Tmb(K) 

CN 

l ->0 

113.491 

10.65 (0.82) 

95.78 (0.37) 

25.788 (1.18) 

0.39 




56.67 (0.59) 

99.31 (0.62) 

92.166 (2.18) 

0.58 




20.75 (0.54) 

156.73 (1.09) 

68.352 (1.92) 

0.28 

H'^CO* 

1^0 

86.754 

2.86 (0.50) 

115.70(10.6) 

111.91 (18.2) 

0.02 

H 

41a 

92.034 

10.75 (0.29) 

103.65 (0.59) 

81.50(1.71) 

0.12 




8.87 (0.39) 

172.86 (7.51) 

352.24 (18.08) 

0.02 

CH3CN 

5/t—>4j.* 

91.987 




<0.015 

N2H* 

1^0 

93.173 

3.85 (0.13) 

93.13 (6.28) 

48.56 (6.28) 

0.07 




1.47 (0.13) 

145.67 (6.28) 

35.74 (6.28) 

0.04 




0.91 (0.13) 

183.60 (6.28) 

30.43 (6.28) 

0.03 

C'^O 

1^0 

109.782 

10.45 (0.60) 

130.29 (3.30) 

111.32 (7.20) 

0.09 

C3H2 

3^2 

145.089 

4.72 (0.20) 

89.20(1.04) 

55.84 (2.74) 

0.08 




2.98 (0.21) 

162.01 (0.31) 

64.91 (4.64) 

0.04 

H2C0* 

2o,2“>1o,i 

145.603 

2.96 (0.25) 

95.22 (0.90) 

37.15 (2.58) 

0.07 

HC3N+H2C0^ 

145.561 

8.20 (0.70) 

93.10(3.30) 

93.54 (10.92) 

0.08 

cs 

3^2 

146.969 

10.96 (2.52) 

96.54 (3.24) 

37.31 (4.49) 

0.27 




7.77 (3.55) 

136.99 (4.27) 

42.03 (16.7) 

0.17 




4.22(1.47) 

182.91 (4.92) 

34.698.89) 

0.11 




W1 (-14" 

-5") 



Line 


Freq(GHz) 

I(Kkms-') 

v,jr(kms *) 

Av(km s“*) 

Tmb(K) 

CN 

1^0 

113.491 

62.38 (0.58) 

103.05 (0.36) 

80.399 (0.87) 

0.73 

H'^CO* 

1^0 

86.754 

3.90 (0.57) 

104.39 (5.51) 

76.888 (13.54) 

0.05 

H 

41a 

92.034 

4.76 (0.18) 

98.51 (1.52) 

80.41 (3.65) 

0.06 

CH3CN 

5 r^ 4 j.* 

91.987 

1.78 (0.17) 

127.93 (3.60) 

71.73 (7.26) 

0.02 

N2H+ 

1^0 

93.173 

6.61 (0.16) 

113.57 (0.74) 

62.23 (1.81) 

0.10 

c'*o 

1^0 

109.782 

24.90 (0.51) 

118.46 (0.65) 

63.32 (1.45) 

0.37 

C3H2 

3^2 

145.089 

13.00 (0.72) 

103.29 (1.55) 

57.80 (3.92) 

0.21 

H2CO 

2o,2“>1o,i 

145.603 

6.12(0.51) 

110.67 (1.18) 

47.71 (2.57) 

0.12 

HC3N 

16^15 

145.561 

3.66(1.52) 

100.41 (18.4) 

93.16(43.1) 

0.04 

CS 

3^2 

145.760 

29.52 (0.26) 

116.01 (0.24) 

54.22 (0.55) 

0.51 


'The PdBI images have been convolved to the same angular resolution (5.9" x 5.9"). ‘The k=0,l,2 and 3 components are blended. The velocity of 
the Gaussian fit is calculated relative to the frequency of the k=0 component, y=91987 MHz. ^The ~326 km s“* component of H 2 CO overlap with 
the ~270 km s“' component of HC 3 N. ^ A possible 180 km s“' component of H 2 CO would overlap with the ~93 km s“' component of HC 3 N. 
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Fig. A.l. Spectral maps of the CN 1—>0 line convolved to an angular resolution of 3.8". To increase the S/N ratio, the velocity reso¬ 
lution has been degraded to 12.8 km s '. The number in the right-upper corner indicates the central channel velocity. Contour levels 
are 0.17 (5xcr) to 1.955 in steps of 0.255 K. Crosses indicate the dynamical center of the galaxy and the position of SN 41.9-1-58. 
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9^66“6B* 60* 

RA (JROOO) 


Fig.A.2. The same as Fig. A.l for the CS 3—>2 line. Contour levels are 0.065 (3xcr) to 1.495 in steps of 0.13 K. The eastern knot at 
velocities ~144-192 km s ' is due to the emission of the H(35 )q' line. 
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Fig. A.3. The same as Fig. A.l for the H(41)a line. First contour and step are 0.028 K (~3xcr). 
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